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We study gravitational perturbations in the Randall-Sundrum two-brane 
■ - - , background with scalar-curvature terms in the action for the branes, allowing 

I for positive as well as negative bulk gravitational constant. In the zero-mode 

■ approximation, we derive the linearized gravitational equations, which have 



(N 



the same form as in the original Randall-Sundrum model but with different 
expressions for the effective physical constants. We develop a generic method 
for finding tachyonic modes in the theory, which, in the model under consid- 
I eration, may exist only if the bulk gravitational constant is negative. In this 

case, if both brane gravitational constants are nonzero, the theory contains 
. one or two tachyonic mass eigenvalues in the gravitational sector. If one of 

I , the brane gravitational constants is set to zero, then either a single tachyonic 

! mass eigenvalue is present or tachyonic modes are totally absent depending 

\ on the relation between the nonzero brane gravitational constant and brane 

separation. In the case of negative bulk gravitational constant, the massive 
gravitational modes have ghost-like character, while the massless gravitational 
mode is not a ghost in the case where tachyons are absent. 
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I. INTRODUCTION 

Linear gravitational perturbations of the flat braneworld models were studied beginning 
from the seminal papers by Randall & Sundrum [1,2], where their spectrum was shown to 
contain, besides the zero mode, also an infinite tower of Kaluza-Klein massive modes. Since 
then, perturbations in various types of braneworld scenarios and in various approximations 
were considered in the fiat case as well as on the cosmological background (which is a much 
more complicated issue still far from being well understood; see [3] and references therein). 
In this paper, we study certain aspects of linear perturbations on a particular simple back- 
ground, which, to our knowledge, have not been previously investigated. We consider the 
Randall-Sundrum two-brane model (the so-called RSI model [1]) supplemented by scalar- 
curvature terms in the action for both branes. Historically, the induced scalar-curvature 
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term for the brane was introduced in [4] (see also [5]) as a method of making gravity on the 
brane effectively four-dimensional even in the flat inflnite bulk space. The corresponding 
cosmological models were initiated in [5,6]. The effects of these scalar-curvature terms in 
the two-brane setup under consideration were recently studied in [7-11]. Perturbations in 
the one-brane (RS2 [2]) counterpart of this model were previously investigated in [12,13]. 

The study of the brancworld models is usually conflncd to the case of positive bulk and 
brane gravitational constants. (The signs of the gravitational constants in this case are 
defined relative to the signs of the conventional matter Lagrangians on the branes.) This 
can be explained by the fact that only positive values of the bulk gravitational constant are 
allowed in the Randall-Sundrum model, since its sign coincides with the sign of the effective 
Newton's constant in that model [1,2]. When one adds curvature terms in the action for 
the branes, this assumption might be relaxed, and, in this paper, we allow for positive as 
well as negative bulk gravitational constant while keeping gravitational constants on the 
branes positive. More generally, one could consider various relations between the signs of 
the four-dimensional gravitational constants in the action for the branes, bulk gravitational 
constant, and conventional matter Lagrangians [8-11]. Although we are not aware of any 
fundamental multidimensional theory that can produce different signs of the gravitational 
couplings for the bulk and branes, we consider this possibility from the viewpoint of the 
effective action regardless of the unknown underlying theory. 

The study of both signs of the bulk gravitational constant is partially motivated by the 
existence of braneworld models with interesting behaviour that require negative brane ten- 
sion. One of them is the model of disappearing dark energy (DDE) [14,15]. The DDE model 
is a braneworld model of expanding universe which, after the current period of acceleration, 
re-enters the matter-dominated regime continuing indefinitely in the future. The merit of 
this model of dark energy is the absence of the cosmological event horizon owing to the fact 
that the universe becomes flat, rather than Dc Sitter, in the asymptotic future. The DDE 
model is based on the generic braneworld action with the bulk cosmological constant and 
brane tension satisfying the Randall-Sundrum constraint, and also including the curvature 
term in the action for the brane. For the consistency of this model with the current cosmo- 
logical observations (specifically, for the condition < 1 on the dark-matter cosmological 
parameter), the brane tension has to be negative [14]. Negative brane tension is required also 
for the existence of unusual 'quiescent' singularities [16] in the AdS-cmbedded braneworld 
models, which occur during the universe expansion and are characterized by finiteness of 
the scale factor, Hubble parameter, and matter density. 

The bulk gravitational constant enters the homogeneous cosmological equations on the 
brane in even power; therefore, its sign does not matter on the level of the homogeneous 
cosmology on the brane [14]. However, the relation between the sign of the bulk gravitational 
constant and the sign of the brane tension is of importance for the small-scale gravitational 
physics in the braneworld, in particular, for the behaviour of cosmological perturbations. 
This can be seen already from the fact that the property of 'locahzation' of five-dimensional 
gravity in the neighbourhood of the brane (that the warp factor locally decreases as one 
moves away from the brane) requires the brane tension and the bulk gravitational constant 
to be of the same sign. This localization property may turn out to be important for a 
consistent braneworld theory and, therefore, for negative-tension branes, it may require 
negative gravitational constant in the bulk. Thus, it seems important to keep open this 
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possibility when generalizing the Randall-Sundrum model by including the curvature terms 
in the action for the branes. 

The issue of ghosts in a theory with positive bulk gravitational constant but arbitrary 
signs of the brane gravitational constants was recently under investigation in [9] , and certain 
regions of parameters were ruled out. In this paper, we consider somewhat complimentary 
situation where the brane gravitational constants are always positive with respect to the 
conventional matter Lagrangian, while the bulk gravitational constant can be of any sign.^ 
We will be mainly concerned with the issue of tachyons in the theory of this kind, which, to 
our knowledge, was not discussed in the literature before. Our original method of finding 
tachyonic modes in the two-brane background is quite general and can be applied to other 
braneworld models leading to the same system of equations for the bulk gravitational modes, 
in particular, to the model with the Gauss-Bonnet action in the bulk considered in [8]. 

It will be shown below that, unlike in the pure Randall-Sundrum case, the presence of 
the curvature terms in the action for the branes leads to a possibility of unwanted tachyonic 
modes if the bulk gravitational constant is negative. We demonstrate that, in this case, 
there can be only one or two tachyonic mass eigenvalues in the theory under consideration 
and determine the range of parameters for which tachyonic modes do not exist. We will also 
see that the massive gravitational modes have ghost-like character in the case of negative 
bulk gravitational constant, while the massless gravitational mode is not a ghost in the case 
where tachyons are absent. These results are in agreement with those of [9]. Thus, it is 
mainly the presence of ghosts in the Kaluza-Klein massive spectrum of gravity that makes 
models with negative bulk gravitational constant problematic. 

This paper is organized in the following manner. After describing the model, we review 
a suitable theory of linear gravitational perturbations on the Randall-Sundrum two-brane 
background. Then, in the zero-mode approximation, we derive the linearized system of 
gravitational equations with matter confined to the visible and/or hidden brane. These 
equations have the same form as in the original RSI model but with different physical con- 
stants. After that, we specially investigate the case of negative bulk gravitational constant 
and show that the linearized theory can contain tachyonic gravitational modes. In this case, 
one or two tachyonic mass eigenvalues are observed if both brane gravitational constants 
are nonzero. If one of the brane gravitational constants is zero, then either a single tachy- 
onic mass eigenvalue is present or tachyonic modes are totally absent; which of these two 
possibility is realized depends upon the values of the nonzero brane gravitational constant 
and brane separation. We determine the range of parameters for which tachyonic modes 
are present or absent. Following [9], we also consider the ghost modes for the radion and 
graviton in our theory. In the case of negative bulk gravitational constant, the range of 
parameters where the radion is ghost free is rather narrow, and all massive gravitational 
modes have ghost-like nature. 

We also calculate the effective gravitational potentials of static matter sources on the 
visible and hidden brane in the generic case. 



The results obtained in one of such theories are usually easy to apply to another theory by using 
the change of the overall sign of the action. 
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II. THE MODEL 



The action of the model in the neighbourhood of one brane has the form 



I (7^-2A)-2/ k\ + I (m^R-2a)+f L{KbA) . (1) 

Vbulk J brane J J brane ^ ' J brane 



where the first part, proportional to the cube of the bulk Planck mass M'^, describes the 
bulk bounded by the brane, and the remaining integrals are taken over the brane. Here, TZ 
is the curvature scalar in the bulk, R is the curvature scalar of the induced metric hah on the 
brane, K is the trace of the tensor of the intrinsic curvature Kah of the brane with respect 
to the inner normal, L {hob, (f>) is the Lagrangian of the matter fields (f) on the brane, and 
integration in (1) implies natural volume elements in the bulk and on the brane. The action 
is similar in the neighbourhood of the other brane. In principle, two branes in our model 
may have different Planck masses m, and, to allow for solutions with flat vacuum branes, 
their tensions a must have opposite signs and satisfy the well-known constraint [1,2] 

Note that, in this paper, we allow for positive as well as negative signs of the bulk Planck 
mass parameter M, the consequences of which will become clear later. 

Action (1) leads to the bulk described by the usual Einstein equation with cosmological 
constant: 

Qab + ^9ab = , (3) 

while the field equation on the brane is 

rri'Gab + (Thab = r^b + (K^b - KbK) , (4) 

where Tab is the stress-energy tensor on the brane stemming from the last term in action 
(!)• 

By contracting the Gauss identity 

Rabc'' = hjhh'hj'h''^nfgk' + KacK," - Kb.Ka" (5) 

on the brane and using Eq. (3), one obtains the 'constraint' equation 

i? - 2A + KabK"-^ (6) 
which, together with (4), implies the following closed scalar equation on the brane: 

{R - 2A) + {m^Gab + ahab - Tab) (m^G'^^ + a/i"' - r"') - ^ (m'R - 4a + r)' = , (7) 

where r ~ hf^^Tgh- 

In the case of a vacuum brane {Tab — 0), Eq. (7) takes the form 

(m^ + ^am^) R + m^ [RabR'^' - \r^^ - ^M^^^ = , (8) 
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where Ars is given by Eq. (2). It should be noted that the second term in Eq. (8) has 
precisely the form of one of the terms in the expression for the conformal anomaly, which 
describes the vacuum polarization at the one-loop level in curved space-time (see, e.g., [17]).^ 
Another useful relation is the Codazzi identity 

Da{K\-h\K)^0, (9) 

which is valid at any timelike hypersurface in the bulk, in particular, on the branes, due 
to Eq. (3). Here, Da denotes the unique covariant derivative on the timelike hypersurface 
associated with the induced metric hat- 

The gravitational equations in the bulk can be integrated by using Gaussian normal 
coordinates, as described, e.g., in [18]. Specifically, in the Gaussian normal coordinates 
{x, y), where x — {x"} are the coordinates on the brane and y is the fifth coordinate in the 
bulk, the metric is written as 

ds^ = dy"^ + hap{x, y)dx'^dx^ . (10) 

Introducing also the tensor of extrinsic curvature Kab of every hypersurface y = const, one 
can obtain the following system of differential equations for the components ha/s and K'^p: 

= R% - KK% - -6% (R + 2A + K^^,K\ - K^) 

^R'^p-KKy-^SyA, (11) 



dhai 
dy 



2ha^K^p, (12) 



where i?"^ are the components of the Ricci tensor of the metric hap induced on the hyper- 
surface y — const, R — R°'a is its scalar curvature, and K — K°'a is the trace of the tensor of 
extrinsic curvature. The second equality in (11) is true by virtue of the 'constraint' equation 
(6). Equations (11) and (12) together with the 'constraint' equation (6) represent the 4+1 
splitting of the Einstein equations in Gaussian normal coordinates. The initial conditions 
for these equations are defined on the brane through Eq. (4). 



III. LINEAR PERTURBATIONS 

Linear perturbations of the RSI model are well studied (see, e.g., [19,20] and references 
therein). Here, we would like to see the modifications arising from the presence of the 
scalar-curvature terms in the action for the branes (nonzero values of the masses m and 



^It is interesting that, while the conformal anomaly term RabR'^'^ — ^R^ cannot be obtained by the 
variation of a local four-dimensional Lagrangian, the very same term is obtained via the variation 
of a local Lagrangian in the five-dimensional braneworld theory under investigation [14] . Also note 
that this term is absent in the original Randall-Sundrum model which has m = 0. 
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m*). Several aspects of this setup were also studied in [8-11]. Our treatment in this and in 
the subsequent section is similar to that of [20]. 

The perturbed metric of our solution in Gaussian normal coordinates has the form (10) 
with 

hap{x,y) ^ (y) [rja/s + la/six, y)\, (13) 

where 

a{y) = exp{-ky) , k = , (14) 

and we emphasize that k can be positive as well as negative depending on the signs of M 
and a. The perturbations of the tensor of extrinsic curvature and of the Einstein tensor 
have the form 

where 7^/3 = 7^/3 — |^a/37; 7 — I'^ai ^-nd □ = d'^da- Here and below, the indices of and 
7q/3 are raised and lowered with respect to the fiat metric ry^g. 

Using the freedom of choice of the coordinates x°' on the brane, one can choose the 
harmonic gauge in which 9"7a/3 = on one of the branes. In this gauge, we have 

Gap = -\^lap (16) 

on that brane. 

In the unperturbed solution, the first (visible) brane is assumed to be at y = 0, and 
the second (hidden) brane is aX y = p. First, we consider the situation where the hidden 
brane does not have matter on it (stress-energy tensor equal to zero). Then, when studying 
perturbations, it is convenient to choose Gaussian normal coordinates with respect to the 
hidden brane. Thus, the hidden brane remains at y = p, while the position of the visible 
brane is linearly perturbed to become y = (f){x), which is the so-called radion degree of 
freedom. Let m and a denote the Planck mass and tension of the visible brane, and let those 
of the hidden brane be m* and cr* = —a, respectively. The linearly perturbed boundary 
equation (4) on the second (hidden) brane becomes 

-mlGy = M'6S^fS: (17) 

where S^/s = K^^ — S^isK, and we have taken into account that the extrinsic curvature 
is calculated with respect to the normal in the positive direction of y. Choosing harmonic 
coordinates on the hidden brane, for which (16) is satisfied, we have 

where a* = a{p) — e'^'' and 7 = 7"a- 

Linearization of the vacuum constraint equation (8) implies the condition 07 = on the 
hidden brane if 
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M^+'^a^ml^O, (19) 

which we assume to be the case. This condition and Eq. (18) imphes the condition d^/dy = 
at the hidden brane. Then the Codazzi relation (9) imphes the condition d {d^^ap) /dy — 
at the same brane. 

Now we turn to Eqs. (11) and (12). Using (15), we can write the second-order differential 
equations for perturbations ja/3 in the bulk. First, we verify that the Gaussian normal 
coordinates remain harmonic in the bulk. We introduce the quantity 

= , (20) 

which is an indicator of the harmonicity of the coordinates on the hypersurface y = const. 
Then we can write the system of differential equations for v°' and 7 that stems from system 
(11), (12): 

^=4A:^ + A;^^, ^ = -J-(2a„^« + n7) + 8fc^, (21) 
oy^ ay ay oy^ oy 

with the following boundary conditions at the hidden brane {y — p): 

„" = 0, ^ = 0, D^^O. gl = 0. (22) 

The unique solution of system (21) in the bulk with the boundary conditions (22) is 

T;"(x,y) = 0, 7(x,y)=7(x), □7(2:) = . (23) 

In particular, this means that the Gaussian normal coordinates which we are using remain 
harmonic with respect to x all over the bulk. 

Taking into account relations (23), from (11), (12) one obtains the system of equations 
for perturbations in the bulk: 

with the boundary condition at y = p which stems from (18): 
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To obtain the boundary equations on the visible brane, one must take into account its 
'bending' in the bulk: y = (f){x). The induced metric on the visible brane in the linear 
approximation becomes 

/l5=(l-2fc0)77„;3 + 7a/3, (26) 

SO that its perturbation is 

75 = 7a/3 - '^kHa(3 , 7$ = 7a/3 + '^kcf>r]af3 ■ (27) 

7 



Substituting it to the boundary condition at the visible brane 

m^G"^ = M^SS'^p + r"^ , (28) 
we obtain the boundary condition at y = 0: 

m^G^p = -^OTa/s + am^/c {d^dp - rj^pU) = r^^ + (77^/3 □ - d^dp) + ■ 

(29) 

Taking trace of this equation, we obtain the equation for the radion field 0: 

-3An(/) = T, (30) 

where r = r]"^Tai3 is the trace of the stress-energy tensor, and A = + 2kw?. Thus, 
the radion field is coupled to the trace of the stress-energy tensor, as is the case in the 
Randall-Sundrum model [19], but with different coupling constant. Using Eq. (30), from 
(29) we obtain 

-^□7,^ = T^p - ^n^pr - Adadpcl> + \m'^ . (31) 

Now we have to solve the bulk equations (24) with the boundary conditions (30), (25) 
and (31). Proceeding to the Fourier transform with momenta pa in the coordinates a;" and 
omitting the tensor indices, we have for the Fourier image ■0(g, y) of 'japi^, u)'- 

il}" - Akil^' + q^e^^^il^ = , (32) 

where the prime denotes the derivative with respect to y, and q — ^—p^ (here we assume 
p2 ^ p'^p^ < 0; the tachyonic case will be studied in Sec. V). After the standard change of 
variable and function 

z{y) = ^ , m = z'xiz) , (33) 



we get the equation 



" + zx'+{z'-4)x = 0, (34) 



in which the prime denotes the derivative with respect to z. Note that z{y) is a monotonic 
function of y for both signs of k, but the sign of z coincides with the sign of the constant k. 
The boundary conditions follow from (31) and (25): 



2T 



-m^ez^oXizo) = ^ + kM^ Uox'izo) + 2x(^o) I , ^0 = ^(0) = ^ , (35) 
Zq k 

mlkhlxiz*) = kM^ [z.x'iz.) + 2xiz.)] , z, = z{p) = ^e'^f , (36) 
where T stands for the Fourier transform of the expression 
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^afs = Tafs - :^r]afsr - Ada 



(37) 



with tensor indices omitted. 

The general solution of the Bessel equation (34) is given by 



xiz)^PJ2{\z\)+QY2 



2 1^1 



(38) 



where J2 and Y2 arc the Bessel functions, P and Q are constants, and the modulus of z 
reflects the fact that the domain of z is positive or negative depending on the sign of k. 
Using the recurrence relations 



zJ!,{z) + 2J2{z) = zJi{z) , zY^{z) + 2Y2{z) = zYi{z) , 
we obtain from (35) and (36): 

2T 



2 ; 2 2 

—m k Zq 



2 1 2 2 

m± z^ 



PJ2{W\)+QY2{W 

PJ2{\zA)+QY2{\z,\) 



+ kM^\z^\ [PJi{\zo\) + QY^[\zo\) 

kM^\z,\ [PJi[\z,\) + QYi{\z,\) 



Zo 



(39) 

(40) 
(41) 



solving which, one finds the constants P and Q and obtains the solution for ip{z): 



QTz^\C*yJ2{\z\) - C:jY2{\z\) 



(42) 



where the constants are given by 



C'y = Y,{\zo\) + ^,k\zo\Y2{\zo\) , C'j = J,{\zo\) + ^A;|zo| ^2(^0!) , 



M3 



CP = Y,{\z.\) - ^,k\zAY2{\z.\) , C} = J,{\z.\) - ^,k\z.\J2{ 



(43) 



These results differ from the similar results [19,20] of the RSI model by the presence of the 
terms containing the brane Planck masses m and in Eqs. (43). 

The spectrum of the model is determined by the equality of the denominator of (42) 
to zero. Introducing the dimensionless variable s — q/\k\ and parameters = km'^/M^, 
fjL^ — kml/M^, and a — e'^f, we obtain the following equation for the spectrum: 



where 



Fi{s) + iJ,sF{s) + aiJ,^sF^{s) + ann^s^ F2{s) = , 



Fi{s) = Ji{s)Y,{as) - Ji{as)Y,{s) , 

F{s) = J2is)Yi{as) - Ji(«s)F2(s) , 
F^s) = J2{as)Yi{s) - Ji{s)Y2{as) , 
^2(5) = J2{as)Y2{s) - J2{s)Y2{as) . 



(44) 



(45) 



If both masses m and are nonzero, then the ultraviolet asymptotics of the spectrum 
for the Kaluza-Klein modes is determined by the zeros of the last term in (44), so that 
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Sn , n > 1 , 46 

a — I 

which coincides with the asymptotics of the spectrum in the Randall-Sundrum model, de- 
termined by the zeros of the first term in (44). 

If m 7^ 0, = 0, then the asymptotics of the spectrum is determined by the second 
term in (44): 

vrn — f , , 

Sn T ' ^ > 1 ■ (47) 

a — i 

If m = 0, m* 7^ 0, then it is determined by the third term in (44): 

vrn + f , , 

Sn f , n > 1 . 48 

a — 1 

In these last two cases, the spectrum is somewhat shifted. 



IV. LINEARIZED EQUATIONS IN THE ZERO-MODE APPROXIMATION 



In the zero-mode approximation [19,20], one considers the hmit as g — 0. In this limit, 
using Eqs. (29), (30), and (42) and expanding all functions of g = \/—p^ in powers of q, 
we obtain the linearized gravity equation in the case of matter present only on the visible 
brane: 



2k 



Ta0 - 



3/1 



p2 



(49) 



where it should be stressed that Gap is the Einstein tensor of the induced metric (26) on 
the brane, and the constants A and B are given by^ 



A^M^ + 2km^ , B^M^ - 2kml . 
The effective Newton's constant Gn is given by the relation 



A - Be-^^p ' 



(50) 



(51) 



and one should note the extra contribution from the radion in (49) , which involves the trace 
of the stress-energy tensor."^ If /c > 0, this contribution is exponentially suppressed for large 
separations between the branes, kp ^ 1. 



^The constant A is the same as in Eq. (30). 

^In Eq. (49) as well as in similar equations of this section, we formally express the radion field 
through the trace of the stress-energy tensor using Eq. (30), similarly to how it is done, e.g., in 
[20]. 
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If matter is present only on the hidden brane, then it still induces curvature on the visible 
brane [19,20]. In our theory, we obtain the result 



G 



a/3 



2k 

Ae^^P - B 



,-2kp 



PaPiS 
p2 



(52) 



where r*^ is the stress-energy tensor on the hidden brane, and r* is its trace. If both branes 
contain matter, then the results on the right-hand sides of (49) and (52) simply add together. 

A few comments are in order about the obtained results. First of all, in the limit of 
zero Planck masses for the branes, m = = 0, we have A = B = M^, and they turn 
to the results previously obtained for the Randall-Sundrum two-brane model [19,20]. The 
presence of two new mass parameters m and extends the freedom of the model. Thus, 
if the constant B turns out to be sufficiently small, then the scalar contribution to the 
right-hand side of (49), proportional to the trace of the stress-energy tensor, may become 
negligibly small. Note, however, that it is not possible to set either the constant A or 
the constant B exactly to zero in our expressions since, in this case, the theory becomes 
singular. This can be seen, e.g., from Eq. (8), in which the first (linear in curvature) term is 
proportional to A [the same property is observed in the general equation (7)]. In particular, 
the nonzero value of the constant B was already assumed in the linearization scheme [see 
Eq. (19)]. The special cases where either A or S is equal to zero must be studied separately. 
Some results in this direction were recently reported in [10], where it was pointed out that 
the linearized theory possesses some additional symmetry in this case. 

To see how this degeneracy arises in some more detail, we turn to the Gauss identity 
(5) again and, following the procedure first employed in [21], contract it once on the brane 
using equations (3) and (4). We obtain the effective equation on the brane that generalizes 
the result of [21] to the presence of the brane curvature term: 



2a 1 

Gab + ARS^hab = ^^^ab + JpjQab ' -J-W^b , 



where Ars is given by (2), 

Qab 



-EEab - EacE'^b + 7^ ( EcdE"^ - -E'^] hab 



(53) 



(54) 



is the quadratic expression with respect to the 'bare' Einstein equation Eab = rn^Gab — Tab 
on the brane, E = h°'''Eab, and Wab = h'^ah^b^cdef^'^'n-'^ is a projection of the bulk Weyl 
tensor Wabcd to the brane. One can see that all the couphngs in (53), including the effective 
cosmological and gravitational constants, are inversely proportional to the constant A, which 
indicates that the theory becomes degenerate in the case A = 0. In the absence of the 
curvature term on the brane (m = 0), we have A — M^, which brings us to the original 
result of [21]. 

Our second remark is that, unlike in the original Randall-Sundrum model (m = m^, = 0), 
in our theory the sign of the constant M is not fixed by the zero-mode approximation: apart 
from the scalar contribution described by the trace of the stress-energy tensor, matching 
with the general-relativity limit fixes only the sign of the overall constant in (49) and (52). 
In particular, for a sufficiently small absolute value of M, namely, \M^/2k\ -C m^, ml , the 
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sign of M does not matter. In the formal limit M — >■ with A; > (hence, k +00), the 
equation for the visible brane (49) turns to the usual linearized Einstein equation. In the 
formal simultaneous limit M — > and cr so that k = a/M^ is fixed, expressions (49) 
and (52) become 



G. 



1 



a/3 



Tq/3 + 



mle '^^P 
2>m? 



Va(3 



PaPp 
p2 



and 



G. 



aP 



m 



!g2fep _j_ 



mt 



,-2kp 



^aP 



VaP 



PaPp 
p2 



(55) 



(56) 



respectively. However, in the following section we will see that the massive gravitational 
modes in the theory with negative value of M have ghost-like nature. 

Finally, we note that our result does not explicitly contain the constant a but contains 
it only in the combination k = a/3M'^. Therefore, for one and the same effective law of 
gravity (49) and (52), the visible brane can have either positive or negative brane tension, 
depending on the sign of M. In particular, the zero-mode graviton is 'localized' around the 
visible brane {k > 0) even if its tension is negative, provided M is also negative. If A; > 0, 
then, in the limit p — > 00, we pass to the one-brane model in Eq. (49), which has the form 
of the corresponding equation of general relativity. 



V. TACHYONIC MODES AND GHOSTS 

In the original Randall-Sundrum model, negative values of the bulk Planck mass M are 
nonphysical because this leads to negative effective Newton's constant. This can be seen by 
setting m = = in Eq. (51), thus having A = B = in it. If M < 0, then the effective 
Newton's constant is negative for any sign of k, which means that the massless graviton 
becomes a ghost. 

The presence of the curvature terms in the action for the brane relaxes the situation 
with the massless gravitational modes and thus relaxes the necessity of dealing only with 
positive values of M. Negative values of M are of interest in view of some of the braneworld 
cosmological models with negative brane tension, in particular, the model of disappearing 
dark energy [14,15], as discussed in the introduction. However, unlike in the pure Randall- 
Sundrum case, the presence of the curvature terms in the action for the branes leads to a 
possibility of unwanted tachyonic modes and ghost-like character of the massive modes in 
the gravitational sector of the theory if M < 0. In this section, we demonstrate that there 
can be only one or two tachyonic mass eigenvalues in the theory under consideration and 
determine the range of parameters for which tachyonic modes do not exist. We also show 
that the massive gravitational modes have ghost-like character in the case of negative M. 

In looking for tachyonic modes, one needs to solve Eq. (32) for — —p^ < 0, i.e., 

V^" - Akil:' - p^e^^^il: = , > , (57) 

where the prime denotes the derivative with respect to y. After the standard change of 
variable and function 
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z{y) = pe'y/k , ^{z) = z\{z) , (58) 

we get the equation for the new function x{z)- 

z\" + zx'-{z' + a)x-0, (59) 

in which the prime denotes the derivative with respect to z. With matter present on the 
visible brane only, the boundary conditions are similar to (35) and (36): 

2T r 1 V 

m\^zlx{zQ) ^^ + kM^ [zox'izo) + 2x(zo)J , zo = ^(0) = | , (60) 



^2 
^0 



-miezlxiz.) = kM^ [z.x'iz.) + 2x{z.)\ , z, = zip) = ^e"" , (61) 

where T is the Fourier transform of expression (37), with tensor indices omitted. 
Solution of (59) is now given by the modified Bessel functions I2 and K2: 

x{z) = Pl2{\z\)+QK,{\z\), (62) 

where P and Q are constants, and the modulus of z again reflects the fact that the domain 
of z is positive or negative depending on the sign of k. The recurrence relations of type (39) 
are valid also for the modified Bessel functions: 

zl!,{z) + 2l2{z) = zh{z) , zK'^iz) + 2K2{z) = -zK^{z) , (63) 

and we can use them in deriving the solution similar to (42): 

where 

C'j = /i(kol) - ^fcko|/2(ko|) , = K,[\z,\) + ^^k\z,\K,[\zo\) 

2 2 



(65) 



At this point, we note that the restriction to the brane at y = and the limit of brane 
separation p — > cxo brings expression (64) to the form obtained for the one-brane case in 
[13]. Our result generalizes it to the case of two branes with arbitrary sign of the bulk 
gravitational constant. 

Tachyonic modes correspond to those values of p for which the denominator of (64) turns 
to zero: 

C;C°: - = . (66) 

Since the transformation p e~^^ followed by A; — —k and m ^ m^, does not change the 
spectrum of the theory, it is sufficient to study only the case /c > 0. Using the dimensionless 
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variable s = pjk and parameters ^l — km^/M^, = kml/M^, and a — e''f, we obtain the 
equation for tachyonic modes: 



where 



E{s) = Di{s) + ijlsD{s) + aijL^sD^{s) + a^n^s^ D^is) = , 



D{s) = l2{s)K^{as) + h{as)K2{s) , 
D^s) = his)K2ias) + hias)K,{s) , 
D2{s) = l2{as)K2{s) - l2{s)K2{as) . 



(67) 



(68) 



Since a > 1 for A: > 0, all functions in (68) are strictly positive for positive s. This 
implies that tachyonic modes are absent in the case fi, fi^ > 0, or, equivalently, M > 0. 

However, tachyonic modes may be present in the opposite case M < 0. It is possible to 
indicate the corresponding range of parameters where tachyonic modes are present or absent. 
To do this, it is convenient to introduce the following auxiliary function of two variables s 
and s: 



E{s, s) = Di{s) + fisD{s) + an^sD^{s) + a^iJL^s^D2{s) 



(69) 



[to be compared with (67)]. By construction, E{s, s) = E{s). 

First, we consider the case where both and /x* are nonzero (in the present case, they 
are then both negative). Then the equation 



E{s, s) = 

gives the two branches of solutions with respect to s: 



(70) 



s±{s) 



2aixix^D2{s) 

and solving the original equation (67) is equivalent to solving one of the equations 

s±(s) 



s . 



(71) 



(72) 



It is easy to verify that the expression under the square root of (71) is strictly positive for 
positive s so that the two solutions s±{s) exist for all s > and are positive. The asymptotic 
behavior of these solutions for small and large s can easily be found: 
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« /i + A** 



(a^ - 1) /Xyu*s ' 



[a 



0, 



(73) 



jji + ajj,^ ± /i — ail 



const , 



oo . 



(74) 
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From these expressions it is clear that the graph of s+(s) crosses the graph of /(s) = s 
at least once for any values of parameters in the range /i, < under consideration. Thus, 
Eq. (72) has a solution, and at least one tachyonic mass eigenvalue is present in this range of 
parameters. Numerical computation indicates that there is exactly one solution connected 
with the branch in this case. 

It is also clear that the graph of definitely crosses the graph of f{s) = s in the case 
where s'_(0) > 1, or 

>1. (75) 



Again, numerical computation indicates that there is only one tachyonic solution connected 
with the branch s_ (s) in this case. They also indicate that tachyonic modes connected with 
the branch s_(s) are absent in the case of the opposite inequality in (75). 

Thus, tachyonic modes exist for all values of parameters in the range /j, ji^ < 0. In the 
case under consideration, M < 0, one can expect tachyonic modes to be absent only if one 
of the brane Planck masses m or is zero.^ We show that this is indeed the case and 
determine the range of masses and brane separations for which tachyonic modes are absent. 

In the case // < 0, //* = 0, the function E(s, s) given by Eq. (69) takes the simple form 

E{s,s) = D,{s) + iJsD{s) , (76) 
and Eq. (70) has a single solution with respect to s: 

= (77) 

^ ^ nD{s) \n\D{s) ^ ^ 

It can be verified that the function Di{s)/D{s) is convex upwards, so that the equation 
s{s) — s has exactly one solution or no solutions in the range s > depending on the value 
of the derivative s'(0). Specifically, a solution exists if s'(0) > 1, and there are no solutions 
in the opposite case s'(0) < 1. Calculating the derivative s'(0), we obtain that exactly one 
tachyonic mass eigenvalue is present in the theory if 

< i (l - a-') , (78) 

and tachyonic modes are absent if the value of the Planck mass m is sufficiently large, 
namely, if 



> 

~ 2 



I (l - . (79) 



In the limit of infinite separation between branes, a = e^'^ — > cxd, the condition of absence of 
tachyonic modes becomes > 1/2, which coincides with the condition A > 0, where A is 



^Tachyonic modes are obviously absent in the Randall-Sundrum model (m = m* = 0) even in 
the case of negative bulk gravitational constant (M < 0), but this model is already excluded as 
resulting in negative effective Newton's constant on the brane (see the beginning of this section). 
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given by (50). Interestingly, this is also the condition of positivity of the effective Newton's 
constant in the zero-mode approximation (49), (52) in the same limit. 

The case ^ = 0, /X* < is analyzed in quite a similar way. Now the function E{s,s) 
given by Eq. (69) takes the form 

E{s, s) = Di{s) + an^sD^s) , (80) 

and Eq. (70) has one solution with respect to s: 

s(s) ^ -^1^ = ^li^ . (81) 
afi^D^{s) a\fi^\D^{s) 

Again, it can be verified that the equation s{s) = s has exactly one solution in the range 
s > if s'(0) > 1, and there are no solutions in the opposite case s'(0) < 1. Calculating 
the derivative s'(0), we obtain that exactly one tachyonic mass eigenvalue is present in the 
theory if 

|//.|<^(a2-l) , (82) 

and tachyonic modes are absent if the value of the Planck mass m* is sufficiently large, 
namely, if 

|/x*|>^(a'-l) . (83) 

We note that our method of finding the range of parameters where tachyonic terms are 
present or absent is not restricted to the model under investigation and can be used whenever 
the equation for tachyonic modes has the form (67), as is the case, e.g., in the theory with 
arbitrary signs of the brane gravitational constants [8-11] and/or with the Gauss-Bonnet 
action in the bulk [8]. 

The issue of ghosts in the gravitational sector of the complementary theory with positive 
value of M but arbitrary signs of the brane gravitational constants was considered in [9], 
and we apply the results obtained therein to our case. First, we start with the radion. The 
radion degree of freedom in our formalism is connected with the possibility of brane bending 
in the bulk. After identifying the physical degrees of freedom for the radion, one can obtain 
the conditions for its ghost-free character in our model using the results of [9]: 

MM >0, (84) 

which we expressed in terms of our parameters /i — km^/M^ and = kml/M^ restricting 
ourselves to the case k > and taking into account that M can be of any sign. Then the 
conditions of absence of both tachyons and radion ghosts in the case M < following from 
(79), (83), and (84) are 

/^* = 0, 1 (l - 6-2*^") < < ^ , and /x = 0, \fx,\>^[e^^P-l) . (85) 
These conditions on the constants of the theory can be seen to be rather restrictive. 



16 



Following [9] , we can also show that the massive gravitational modes in the theory under 
consideration have ghost-like nature. For free metric perturbations in the form (13), (14) 
described by the transverse traceless modes japix, y) with the boundary conditions (25) and 
(31) in which we set = and Tap = 0, one obtains the gravitational part of action (1) to 
quadratic order in the form 

]\/f CP r 

+ T / ^^^""^^-/^ + ^e-'''^ I dx^'n^ap . (86) 
y=0 y=p 

Expanding the perturbation in the modes '^{q,y) that are solutions of Eq. (32) with the 
corresponding boundary conditions, 

lapix, y)^Yl Xap(q, x)'^{q, y) , (87) 
substituting this expansion into action (86), and using the orthogonality condition 

/' dye-^>'yiP{qi, y)iP{q2, y) + mV(gi, 0)^(^2, 0) + mle-^'^P'il^iqr, p)V(g2, p) = (88) 
for qi q2, we arrive at the following quadratic effective action (cf. with [9]): 

S-\Y.cJ dxx'-^iq, x) (□ - q') Xap{q, x) , (89) 

where 

C, = r dye-^'y[^l^{q, y)f + m2[V^(g, 0)]^ + mle~'''^[il^{q, p)f . (90) 

•J 

For the massless mode (g = 0), we have %l]{Q,y) = const, and the constant Cq is given by 

Co = ^[V'(0, 0)]^ [l + 2/. - (1 - 2//,) e-"'''] (91) 

and is positive in all cases in which tachyonic modes are absent in the theory, as can be seen 
from conditions (79) and (83). Thus, the massless graviton is not a ghost. However, for the 
massive modes (g 7^ 0), using Eq. (32), one can obtain the expression 

C, = — / dye-'''yW{q,y)]\ q^O, (92) 

which shows that the massive modes have positive norm in the case M > 0, and have 
ghost-like nature in the case M < 0. 



VI. CORRECTIONS TO NEWTON'S LAW 

In this section, we compute the gravitational potential V{r) on the visible brane induced 
by a static point source located on the visible or hidden brane and determine corrections to 
the Newton's law in the physically reasonable case A; > 0, i.e., where the zero-mode graviton 
is localized around the visible brane. Our starting formula is 

/^(r^ = -(l + 2y), (93) 
where h^^g^ is the induced metric on the visible brane. 
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A. Matter source on the visible brane 



If matter source is on the visible brane only, then the induced metric /i^™'' is given by 
Eq. (26). The stress-energy tensor of a static point source of mass M. is Too = A4S (f) 
with other components being zero, and r = — Tqo- Its Fourier image is Too(Pa) — 27rA45 (po), 
containing only tachyonic modes. Hence, we use the formulas of Sec. V. 

We take into account (30) and (64) for the Fourier transform of the solution for the in- 
duced metric on the first brane with the source on the same brane. Collecting all expressions 
together, we obtain the Fourier image of the gravitational potential 



' ' f(W 



(94) 



where the function f{p) denotes the second fraction in expression (64) taken at the position 
of the visible brane {z — Zq): 

The potential V{r) is obtained by taking the inverse Fourier transform of (94): 



V{r) ■■ 



where 



I[r) 



AA 



(96) 



/(r) = / dsaiYiikrs) ^(s) , ^(s) = ^ ' 1^* ^ , (97) 

J -C/ ( s ) 

^ ^ 

and E{s) is given by Eq. (67). 

The integral in (97) cannot be evaluated exactly, but it can be approximated in different 
regions of r, as it is done, e.g., in [12,22]. 

1. On very large spatial scales kr ^ a = e^'', we need the asymptotics of the function 
^(s) in the region as -C krs ~ 1. It is given by the expression 

^ ^ 1 + 2/x + a-2 (2/i, - 1) s 

{a^ - 1) (2^, - 1) [3a^ - 1 + 2//, (1 - a^)] + Aa^ In a 
^ 4[a2(i + 2//) + 2//,-l]' ^ ^ 

and the integral in (97) is approximated by using the regularization [12,22] 

oo oo 

j s sm{krs)ds — > ^^^j ssm{krs)e~^^ds — (99) 



with the result 
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7(r) 



TT 



l + 2/x + a-2 - 1) 
The potential in this region has Newtonian form 

GM 



V{r) = 



G = Gn 1 + 



const . 



B \ 



3Aa^ 



) 



(100) 



;ioi) 



where given by (51), which is in complete agreement with the zero- mode approximation 
(49). The theory has continuous Newtonian hmit as a — * oo. 

2. On the scales 1 <^ kr <^ a, it is the region of integration 1/a <^ s <^ 1 that 
substantially contributes to the integral in (97). In this region, we have 



1 1 log(s/2) 



1 + 2/x s (l + 2//)2 
substituting which to (97), we obtain 



/(r) 



TT 



+ 



TT 



and 



V{r) 



GM 



1 + 



1 + 2/x 2(1 + 2/^)2 (A;r)2 
2M3 \ 



3A{kr) 



G = Gn 1 - 



B 



(102) 



(103) 



(104) 



We observe corrections to the Newtonian potential similar to those of the Randall- 
Sundrum model [2] but with somewhat different relative constant. 

3. In the case fcr <^ 1, we can use the asymptotics for the function ^'(s) at infinity: 



IjLS+ 1 + 15/^/8 



(105) 



This case is further partitioned into two asymptotic regions, depending on the magnitude 
of ix. 

3a. kr <C \iJ>\. Here, if /j, is negative, then it cannot be small by absolute value since, in 
this latter case, the theory contains tachyons. We obtain 



7(r) 



n 



+ 



1 /15 



+ 



2/i A* \ 8 jj, 



kr log 



and 

V{r) 



GM kM fib 1\, 
Sr + - log 



37r2m2 I 8 fi 



— H — ] kr 

8 . 



G 



1 A* + 2/3 



87rm2 + 1/2 



(106) 



• (107) 



The logarithmic corrections in these expressions assume that the expression 15/8 + l//xis 
not very small by absolute value. 

3b. II <^ kr <^ 1. Here we must consider only positive /i. In this case. 



'^'^ " Vr ' 



(108) 
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and the gravitational law is five-dimensional: 



The expressions for the gravitational potentials obtained in cases 3a and 3b are analogous 
to those of [4]. The same results would be obtained in the linear approximation in the one- 
brane case considered in [13], as is clear from the remark made after Eq. (65), which identifies 
the corresponding propagators. 



B. Mater source on the hidden brane 

In a similar way one can consider the case where the stationary matter resides on the 
hidden brane with mass A4* defined as Tqq = A4*S{f). In this case, we obtain the following 
expression for the gravitational potential on the visible brane: 



V{r) 



I{r), 



where 



(110) 



/(r) = j dssm{krs)'^{s) , *(s) 



E(s) 



(111) 



and E{s) is given by Eq. (67). 

The asymptotic expressions for ^(s) can be found in various regions: 



2a 



- 1 + a2(2y^+ 1) 



- • - + 0{s) , s^a-\ 



(112) 



27re" 



(1 + 2jji)jji^^/as 



\ 8 n^l as 



, q;-^<s<1, /i*7^0, (113) 



V2 



Txas e 



, q;^<s<1, //* = 0, 



(114) 



2^/ae 



S > 1 , /X^ + ^ , 



HH^as^ + (a* + A**« + 15/x/x*a/8)s 



(115) 



*(s) !^ 2v«e-"% s>l, i_i^i_i^^O. (116) 

Using these expression, it is not difficult to obtain the estimates for the gravitational 
potential V{r) caused by the presence of the static source on the hidden brane in various 
regions. We have 
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V{r) 



kr a , 



(117) 



where Gn is given by (51). Again, in this distance range, the result can be obtained by 
using the zero-mode approximation (52). 

In the case kr <C a, the result crucially depends on whether /j,^ is zero or not: 



V{r) 



cik'^M* 



1-C2 



kr" 



a 



kr <^a, /Li* = , 



(118) 



V{r) 



2k^M, 




C5 + ---)■ I - 




, kr <^a , //* 7^ , 

(119) 



Here the constants c„ take the following approximate values: 

ci«1.77, C2«1.02, C3^l.3, C4«0.35, cg 0.06 , C6«1.12, C7«0.24. 

(120) 

Expression (119) is not valid for sufficiently small //* since wc know that, in the limit of 
/i* 0, the asymptotics changes to (118). In fact, comparison with the exact numerical 
integration of (111) shows that our approximate result (119) is only good for |/Li*| > 0.5. 



VII. DISCUSSION 

It is known that the braneworld model becomes rather rich in its cosmological mani- 
festations if curvature term is present in the action for the brane (see [4-6,12,14-16]). In 
this paper, we studied the model with scalar-curvature terms for the branes on the original 
Randall-Sundrum two-branc background. The linearized gravitational equations (49) and 
(52) in this case have the same structure as in the original Randall-Sundrum model but 
with different physical constants. In the limit of vanishing brane Planck masses m and m*, 
they tend to the known results [12,20,22], which are physically reasonable only for the bulk 
Planck mass M > 0. In the opposite limit of M — > while k = a/M^ is fixed, they produce 
reasonable results (55) and (56) independently of the sign of M. 

In this paper, we developed a general method for detecting tachyonic modes in the 
braneworld theory, which can be generalized to theories giving rise to equations of the 
similar kind, e.g., the theory with Gauss-Bonnet term in the bulk action [8]. In our case, 
for negative values of M, the linearized theory typically contains tachyonic modes in the 
gravitational sector. If both brane Planck masses are nonzero, then we have one or two 
tachyonic mass eigenvalues depending on the constants of the theory, the conditions of 
which were determined in Sec. V. However, in the case where one of the brane Planck 
masses is zero, tachyonic modes are absent if the other brane has sufficiently high Planck 
mass [given by Eq. (79) for the visible brane, and by Eq. (83) for the hidden brane]. In 
the case of negative bulk Planck mass M, the zero-mode graviton is 'localized' around the 
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brane with negative tension and is not a ghost in ah cases where tachyonic modes are absent 
in the theory. However, in all cases with negative M, the massive gravitational modes of 
the theory under consideration have ghost-like character. The conditions of absence of both 
radion ghosts and gravitational tachyons are expressed by (85). 

Exploring both signs of the bulk gravitational constant may be interesting in connection 
with some braneworld cosmological models requiring negative brane tension, such as the 
model of disappearing dark energy (DDE) recently discussed in [14,15] or the braneworld 
models with 'quiescent' cosmological singularities during expansion [16]. The DDE model 
[14,15] represents a cosmological braneworld with the Randall-Sundrum constraint (2), neg- 
ative brane tension, and the condition |//| > 1/2, which is required for physical consis- 
tency and which implies inequality (79). The 'quiescent' singularities in the AdS-embedded 
braneworld models occur during the universe expansion and are characterized by finiteness 
of the scale factor, Hubble parameter, and matter density. The braneworld cosmological 
equations involve the bulk gravitational constant only in even power; therefore, their be- 
haviour is independent of its sign on the homogeneous and isotropic level [14] . However, as 
noted in the introduction, the sign of the bulk gravitational constant is of importance for the 
small-scale gravitational physics in a braneworld universe, in particular, for the behaviour 
of cosmological perturbations. The results of the present paper indicate that models with 
negative bulk gravitational constant can be free from tachyons, although they are plagued 
with massive ghosts in the gravitational sector. Perhaps, the unwanted situation with ghosts 
can be remedied by modifications of the bulk action. It should be emphasized that these 
results do not relate to the braneworld cosmological models with negative brane tension but 
positive bulk and brane gravitational constants, which require future investigation. 
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